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Abstract. We introduce and study a new interacting particles model with 
a wall and two kinds of interactions - blocking and pushing - which maintain 
particles in a certain order. Wo show that it involves a random matrix model. 



Resume. Nous introduisons et etudions un nouveau modele de particules en 
interaction avec mur. Les particules de ce modele se bloquent et se poussent 
les une les autres, afin de rester dans un ordre donne. Nous montrons que ce 
modele est lie a un modele de matrices alatoires. 



1. Interacting particles model 

Let us consider k ordered particles evolving in discrete time on the positive real 
line with interactions that maintain their orderings. The particles are labeled in 
increasing order from 1 to /c. Thus for t S N, we have 

0<Xiit)<---<Xk{t), 

where Xi{t) is the position of the i*^ particle at time t > 0. Particles are initially 
all at 0. The particles jump at times n — and n, n G N*. Let us consider two 
independent families 

{^{hn- ^))t=i.....k-n>i, and "•))i=l,...,fc; n>l, 

of independent random variables having an exponential law of mean 1. For conve- 
nience, we suppose that there is a static particle which always stays at 0. We call 
it the 0*^ particle, and denote Xo{t) its position at time t >0. 
At time n — 1/2, for i = 1, . . . , fc, in that order, the i*'' particle tries to jump to the 
left according to a jump of size ^{i,n— i) being blocked by the old position of the 
{i — 1)'^ particle. In other words : 

• If Xi{n — 1) — ^(1, 71 — i) < 0, then the I''* particle is blocked by 0, else it 
jumps at Xi{n — 1) — ^(1, n — i), i.e. 

Xi{72 - i) = max(0,Xi(n - 1) - ^{l,n - ^)). 

• For i = 2, . . . , fc, if Xi{n - 1) - i{i,n - |) < X,_i{n - 1), then the i*'' 
particle is blocked by Xi^i{n — 1), else it jumps at Xi{n — 1) — £_{i,n~- i), 
i.e. 

Xi{n - ^) = max(Xi_i(n - l),Xi{n - 1) - ^{i,n- i)). 
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At time particles jump successively to the right according to an exponentially 
distributed jump of mean 1, pushing all the particles to their right. The order 
in which the particles jump is given by their labels. Thus for i = l,...,k, if 
Xi_i(n) > Xi{n — 5) then the i*^ particle is pushed before jumping, else it jumps 
from Xi{n — \) to Xi{n — \) + ^{i,n): 

Xi{n) = max{Xi_i{n),X,{n- i)) +^{i,n). 



Xi{n-1) X2(,n-1) Xsin-l) 




Figure 1. An exemple of blocking and pushing interactions be- 
tween times n — 1 and n for fc = 3. 

There is another description of the dynamic of this model which is equivalent 
to the previous. At each time n G N*, each particle successively attempts to jump 
first to the left then to the right, according to independent exponentially distributed 
jumps of mean 1. The order in which the particles jump is given by their labels. At 
time n G N*, the 1** particle jumps to the left being blocked by then immediately 
to the right, pushing the i^^ particles, i = 2, . . . , fc. Then the second particle jumps 
to the left, being blocked by max(Xi(n — 1), Xi{n)), then to the right, pushing the 
i*^ particles, i = 3, . . . ,k, and so on. In other words, for n G N*, i = 1, . . . , fc, 

Xi{n) = max(Xj_i(n),Xi_i(n - l),Xi{n - 1) - ^~{i,n)) +^+(i,n), 

where {i,'n))i^i^,,,^k;nen, and n))i=i,...,fe;„gN are two independent families 

of independent random variables having an exponential law of mean 1. 

2. Results 

Let us denote A4k,m{^) the real vector space oikxm real matrices. We put on 
it the Euclidean structure defined by the scalar product, 

(M, N) = tv{MN*), M,Ne A^fc,„(M). 

Our choice of the Euclidean structure above defines a notion of standard Gaussian 
variable on 7W„,m(IR). We write Ah for the set {M G Mk,k{^) : M + M* = 0} of 
antisymmetric k x k real matrices, and iAk for the set {iM : M G Ak}- Since a 
matrix in iAk is Hermitian, it has real eigenvalues Ai > A2 > • • • > Afc. Morevcr, 
antisymmetry implies that Xk-i+i = — Ai, for i = I,-- - ,[k/2] + 1, in particular 
\k/2]+i = when k is odd. 
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Our main result is that the positions of the particles of our interacting particles 
model can be interpreted as eigenvalues of a random walk on the set of matrices 
iAk+i- 

Theorem 2.1. Let k be a positive integer and (M (n), n > 0), be a discrete process 
on iAk+i defined by 

1=1 ^ ^ 

where the Yi's are independent standard Gaussian variables in A4k+i,2i^)- For 
n€zN, let Ai(ri) be the largest eigenvalue of M{n). Then the processes 

{Ai{n), n > 0) and {Xk{n), n > 0), 

have the same distribution. 

For a matrix M G iAk+i and m£{l,...,fc + l}, the main minor of order m of 
M is the submatrix 

{My , 1 < i,i < m}. 

Theorem 2.2. Let M{n),n > 0, be a discrete process on iAk+i, defined as in 
theorem \2.1\ For n G N, m = 2, . . . + we denote A^™''(n), the largest eigenvalue 
of the main minor of order m of M{n). Then for each fixed n G N*, the random 
vectors 

(Af>(n),...,Ai'=+^^(n)) and . . . , Xfc(n)), 

have the same distribution. 

Let us notice that there already exists a version of theorem 12.11 with particles 
jumping by one (see section 2.3 of [11] ) and a continuous version which involves 
reflected Brownian motions with a wall and Brownian motions conditioned to never 
collide with each other or the wall (see [T]). 

There is a variant of our model with no wall and no left-jumps which has been 
extensively studied (see for instance Johansson [5 , Dicker and Warren [5] , or War- 
ren and Windridge |11)). It involves random matrices from the Laguerre Unitary 
Ensemble. Indeed, in that case, the position of the rightmost particle has the same 
law as the largest eigenvalue of the process (M(n,), n > 0) defined by 

n 

M{n)=^ZiZ'^, n>0, 
1=1 

where the \/2^i's are independent standard Gaussian variables on C*^. Let us 
mention that this model without a wall is equivalent to a maybe more famous one 
called the TASEP (totally asymmetric simple exclusion process) . In this last model, 
we consider infinitely many ordered particles evolving on Z as follows. Initially 
there is one and only one particle on each point of Z_ . Particles are equipped with 
independent Poisson clocks of intensity 1 and each of them jumps by one to the 
right only if its clock rings and the point just to its right is empty. Particles are 
labeled by N* from the right to the left. Then, the time needed for the fc*'' particles 
to make n jumps is exactly the position of the rightmost particle at time n in the 
model without a wall and exponential right jumps. 
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3. CONSEQUENCES 

Thanks to the previous resuhs we can deduce some properties of the interacting 
particles from the known properties concerning the matrices M(n),n > 0. For 
instance, the next proposition fohows immediately from theorem 12.11 and theorem 
5.2 of 3 . 

Proposition 3.1. For n > [^^^], the distribution function of Xk{n) is given by 
P{Xk{n) <t)= det ( /"*a;2j+.+„-[i^]-3+i{. „„,„} e"" dx) j^j^j^j , t e R+. 



Moreover, as n goes to infinity, the process {^M{[nt]),t > 0) converges in 
distribution to a Brownian motion on iAk+i- Since our interacting particles model 
converges to a model of reflecting Brownian motions, proposition 2 of 1^ follows 
from theorem 12. II 



4. PROOFS 



Theorem [271] is proved by proposition l4.7l Theorem [2T2] follows from propositions 
14.81 and 14.91 Let us denote by </> the function from R to M, defined by (j){x) — 
^exp{—\x\), X £R, and consider the random walk S'„, n > 1, on M, starting from 0, 
whose increments have a density equal to (j). The next three lemmas are elementary. 

Lemma 4.1. The process {\Sn\)n>i is a Markov process, with a transition density 
q given by 

q{x,y)^(f>{x + y) + (f>{x-y), x,yeR+. 

Proof. This is a simple computation which holds for any symmetric random walk. 
□ 

Lemma 4.2. Let us consider a process M{n),n > 0, on iA2, defined by 

1=1 ^ ^ 

where the Yi 's are independent standard Gaussian variables in A^2.2(R)- Then the 
process of the only positive eigenvalue of (Af (n), n > 0) is Markovian with transition 
density given by q. 

Proof. The Yi 's are 2x2 independent real random matrices whose entries are inde- 

/y Y 

pendent standard Gaussian random variables on M. Let us write Yi = 



The matrix Yiuj{l)Yi* is equal to 

liYi^i^i - Yi^2Yi.3) 

^^{Yi,2Yi,3 - niYi^i) 

For a S M, we have 

]g(^g-iQ(Y,,4yi,i-Y,.2yi.3)^ = ^ 



Thus, the random variables Yi nYi^i — 11,2^,3) I ~ 1, . ■ ■ ,n, are independent, with a 
density equal to (f>. We conclude using lemma HTT] □ 
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Lemma 4.3. Let r be a real number. Let us consider (^^)n>i and (^,^)n>i two 
independent families of independent random variables having an exponential law of 
mean 1. The Markov process Z{n),n > 1, defined by 

Z{n) = max(Z(n - 1) - C-^) +^rt > 1. 

has transition density 

Pr(x,y) = (f>{x-y) + e'^''(l){x + y), x,y>r. 

Proof. This is a simple computation. □ 

Let us notice that when r = 0, the Markov process {Z{n),n > 0), describes 
the evolution of the first particle. As its transition kernel po is the same as the 
transition kernel q defined in lemma 14. fl theorem 12.11 follows when k = 1 from 
lemma [L2] The general case is more complicated. 

Definition 4.4. We define the function dk : — > R &?/ 

• when k = 2p, p G N* , 

p 

dk{x)=c:^^ n (2^4 

l<i<j<p 4=1 

• when k = 2p - 1, p e W , 

^'^^"^'l c,^ni<.<,<p(-f-^|) else, 

where 

l<i<j<p l<i<p 

The next proposition gives the transition density of the process of eigenvalues of 
the process M(n), n > 0. For the computation, we need a generalized Cauchy-Binet 
identity (see for instance Johansson [6]). Let {E,B,m) be a measure space, and let 

and V'i, 1 !i i^j 1^ "-i be measurable functions such that the (j^iipj^s are integrable. 
Then 

(1) Aet( I (j)i[x)iij{x)dm{x)\ ^ — I det det TT dm(a;fc). 

We also need an identity which expresses interlacing conditions with the help of a 
determinant (see Warren [9J). For x,y € M" we write a; ^ y if x and y are interlaced, 
i.e. 

Xl > yi > X2 > ■ ■ ■ > Xn > yn 

and we write x )~ y when 

Xl > yi > X2 > ■ ■ ■ > Xn > yn- 

When X £ R"+^ and y S M" we add the relation i/„ > a;„+i (resp. j/„ > Xn+i)- Let 
X and y be two vectors in M" such that xi > ■ ■ ■ > Xn and yi > ■ ■ ■ > j/n- Then 

(2) ^xyy — det(l{a;.>y^})„xn- 
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For k> 1, we denote by the subset of defined by 

Ck = {x G RIt^I : xi> ■ ■ ■ > a;[fc+i] > 0}. 

Proposition 4.5. Let us consider the process {M{n), n > 0), defined as in theorem 
\2.1[ For n e N, let A{n) be the first [-^y^] largest eigenvalues of M{n), ordered such 
that 

Ai(n) > • • • > \h+l]{n) > 0. 
Then {A{n),n > 0), is a Markov process with a transition density Pk given by 

Pk{\P) = ^det(</)(A,-/3,) + (-l)'=+V(A.+/?,))i<„.<.i±i,, A,/3eCfe. 

Proof. The Markov property follows from the fact that the matrices 

^'(-. o)^*' 

are independent and have an invariant distribution for the action of the orthogonal 
group by conjugacy. Proposition 4.8 of [3 ensures that the transition density with 
respect to the Lebesgue measure of the positive eigenvalues of M{n), n > 1, is given 

by 

• when k = 2p, 

• when fc = 2p — 1 , 
where 

When k is even, using identity ([2]), we write l{A,,3;^z}e~ ^^=i'''''^^'~^^''' as 

det(l^,<A,e~(^^-^'))pxpdet(l,.<^^e^(ft-^-))pxp, 
and use the Cauchy-Binet identity to get the proposition. 

When k is odd, we introduce the measure fi on R, defined hy fi ~ 6q + m, where 
So is the Dirac measure at and m is the Lebesgue measure on R. We have the 
identity 

(l){x - y) + (l,{x + y) ^ / l[o,.^j,]e-("+^-2.)^^(^)^ 
Jr 

Thus using the Cauchy-Binet identity with the measure we get that the determi- 
nant of the proposition is equal to 

4 / det(l[o,A,](z.)e-^^^-^-0det(l[o,;3,](^.)e-(^^-^'0 f[ dfi{z,^). 

Using identity ([2]), we obtain that it is equal to 

•^K" m=l 
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We integrate over Zp in the last integral and use the fact that the coordinates of j5 
and A are strictly positive to get the proposition. □ 

The next proposition gives the transition density of the Markov process n > 

0) defined in section [T] For m > 0, we denote (Z)*^™^ the m*^ derivative of 0, and we 
define the function by 

r+oo 1 



We easily obtain that 



,(")(a;) = J 2 



i(-l)"e-^ ifa;>0 



^ else. 



1 -Eiy ^""^''"'^ else. 



For fc > 2, we denote by Vk the subset of M.'^ defined by 

Vk = {x e R'' : < xi < X2 < ■ ■ ■ < xk}. 

Proposition 4.6. The Markov process X{n) = {Xi{n), . . . , Xk{n)),n > 0, has a 
transition density Qk given by 

Qk{y,y ) = dei{ai^j{ijr,y'j))i<ij<k, y,y G Vk, 
where for x, x' € R, 

Proof. Let us show the proposition by induction on k. For fc = 1, the equality 
holds by lemma 14.21 Suppose that it is true for fc — 1. We write Ci, . . . , Cfc, and 
Li, . . . , Lk for the columns and the rows of the matrix of which we compute the 
determinant. There are two cases : 

• If yk_i > yfe-i then for i = l,...,fc— l,yi< y'k^i < y'k- ^ quick calculation 
shows that all the components of the column Ck + e^''-^^'^''Ck-i are equal 
to zero except the last one, which is equal to Py' iVk-i V ykiy'k)^ where Pr 
is defined in lemma l473l for r S R. 

• If y'k-i < 2/fc-i then for i = 1, . . . , fc - 1, < yk^i < yk- We replace the 
last line Lk by the line Lk — e^''^'-^^'' Lk-i having all its components equal 
to zero except the last one, which is equal to Py^^-^ {yk, y'k)- 

Then we conclude developing the determinant according to its last column in the 
first case or to its last row in the second one, and using the induction property. □ 

We have now all the ingredients needed to prove theorems 12.11 and 12.21 For A 
an element of the adherence Ck of Cfc, we denote GTfc(A) the subset of R[^~1[t~1 
defined by 

GTfc(A) = {(a;(2),... ,x^''+^^) : x^''+^^ = X,x'^'^ G m}^\ h x^'-^\3 < i<k+l}. 
We let 

GTfe=U;,gc.GTfe(A). 
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If (a;(2),...^a;('=+i)) is an element of GT^, then \ • • • ,a;f+^^_) belongs to the 
adherence Vk of Vk- Thus we define as the Markov kernel on Ck x Vk such that 
for A € Ck, the probability measure Lfe(A, .) is the image of the uniform probability 
measure on GTfc(A) by the projection p : GTk — > 2?^ defined by 

where ,a;('=+i)) G GT^. For A G Cfe, the volume of GTfc(A) is given by 

dfc(A). Thus Lj,(A, .) has a density with respect to the Lebesgue measure on GTj, 
given by 

Rogers and Pitman proved in |5] that it is sufficient to show that the intertwining Q 
holds, to get the equality in law of the processes (Ai(n), n > 0) and {Xk{n), n> 0). 
So theorem 12.11 follows from proposition 14.71 



Proposition 4.7. 

(3) LkQk = PkLk 

Proof. The proof is the same as the one of proposition 6 in [1]. We use the de- 
terminantal expressions for Qk and Pk to show that both sides of equality ^ are 
equal to the same determinant. For this we use that the coefiicients a^j 's given in 
proposition 14.61 satisfv for x,x' € ]R-|_, 

aij{x,x') = / ai-ij{u, x') du, 

r+oo 

ai,j{x,x) = - aij+i{u,x')du, 

J x' 

a2^2j (a:^, 0) = 0, 02^,2*-! (0, x) = 1, 02*^ (0, x) = 0, 2i < j 

The computation of the left hand side of ([3]) rests on the first identity. The com- 
putation of the right hand side rests on the others. □ 

Let us notice that there are many reasons to think that we could construct an 
enlarged process ((F*^^^ (n), . . . , F^'^"'"^^ (n)), n > O), living on GT^, such that the 

processes {J^^^^^{n),n > 0), and ((r/^^(n), . . . , r/''+^^(n)), n > O), would have 
respectively the same law as {A{n),n > 0) and {Xi(n), . . . , Xk(n),n > 0). This 
would imply theorem 12. II 

The measure Lfe(0, .) is the Dirac measure at the null vector of GTfe(O). Thus, 
the following proposition is an immediate consequence of proposition 14.71 

Proposition 4.8. 

Q^(0,.) = P,-L,(0,.) 

Keeping the same notations as in theorem l2.2[ we have the following proposition, 
from which theorem 12^ follows. 



Proposition 4.9. PJ^Lk{0, ) is the law of the random variable 

{Kf\n),...,K^t+'\n))). 

Proof. The density of the positive eigenvalues A(n) of M{n) is given by Pk{0, ■ 
Then the proposition follows immediately from theorem 3.4 of [3]. □ 
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5. CONCLUDING REMARKS 

As recalled, the model with a wall that we have introduced is a variant of another 
one with no wall and no left-jumps. For this last model, the proofs of the analogue 
results as those of theorems 12.11 and 12.21 rest on some combinatorial properties 
of Young tableaux. Indeed, Young tableaux are used to describe the irreducible 
representations of the unitary group. The matrices from the Laguerre Unitary 
Ensemble belong to the set of Hermitian matrices which is, up to a multiplication 
by the complex i, the Lie algebra of the Unitary group. Their laws are invariant 
for the action of the unitary group by conjugacy. It is a general result that the 
law of their eigenvalues can be deduced from some combinatorial properties of the 
irreducible representations of the unitary group. 

In our case, the distribution of the eigenvalues of (A/(n), n > 0) can be deduced 
from combinatorial properties of the irreducible representations of the orthogonal 
group (see ^ for details). Many combinatorial approaches have been developed 
to describe these representations. Among them we can mention the orthogonal 
tableaux and the analogue of the Robinson Schensted algorithm for the orthogonal 
group (see Sundaram jlOj). or more recently those based on the very general theory 
of crystal graphs of Kashiwara [7] . None of them seems to lead to the interacting 
particles model with a wall that we have introduced. It would be interesting to 
find what kind of tableau involved in the description of irreducible representations 
of the orthogonal group would lead to this model. 
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